The chemical bond in external electric fields: energies, geometries, and vibrational Stark shifts of diatomic molecules.

It is shown that the response of molecular properties of diatomics such as the total energy, the bond length, and the vibrational Stark shift to an external homogenous electric field (EF) can be predicted from field-free observable properties such as the equilibrium bond length, the bond dissociation energy, the polarizability and dipole moment functions, and the vibrational frequency. Delley [J. Mol. Struct.: THEOCHEM 434, 229 (1998)] suggested to approximate the potential energy surface under an EF by a Morse function augmented with a EF term proportional to the internuclear separation. In this work, this term is replaced by the expression of the field-induced energy change which yields a field-perturbed Morse potential that tends to a constant asymptotic limit when the EF term itself become proportional to the sum of the polarizabilities of the separated atoms. The model is validated by comparison with direct calculations on nine diatomics, five homo-nuclear (H2, N2, O2, F2, and Cl2) and four hetero-nuclear (HF, HCl, CO, and NO), covering a range and combinations of dipole moments and polarizabilities. Calculations were conducted at the quadratic configuration interaction with single and double excitations (QCISD) and density functional theory (DFT)-B3LYP levels of theory using the 6-311++G(3df,2pd) basis set. All results agree closely at the two levels of theory except for the Stark effect of NO which is not correctly predicted by QCISD calculations as further calculations, including at the coupled cluster with single and double excitation (CCSD) level of theory, demonstrate.